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The need for an accurate and realistic description of inelastic, high- 
temperature behavior of materials has attracted the attention of numerous 
investigators in recent years. As a result of their concerted efforts, a number of 
constitutive models called viscoplastic models has emerged. These models consider 
all the inelastic strain (including creep, relaxation, etc.) as a unified quantity 
and, in general, provide a better description of high-temperature inelastic behavior 
of materials. Since these models incorporate as much material science as possible, 
their mathematical structure is very complex. To use these models effectively and 
reliably for the analysis and design of turbine engine hot section components, we 
must demonstrate their feasibility by incorporating the models in nonlinear finite 
element structural analysis codes and perform nonlinear stress/life analyses for 
these components. 

With this objective in mind, the NASA Lewis Research Center has mounted an 
in-house effort to implement some of the more commonly used viscoplastic models into 
the general purpose finite element structural analysis code - MARC. The intent is to 
provide the designer and structural analyst with a tool for realistic and rational 
analysis and designs of hot section components. 

This paper gives a brief description of the implementation in MARC of two 
viscoplastic models developed by Robinson (refs. 1 and 2). One model is for 
isotropic materials and the other for metal matrix composites. Also presented are 
analytical results obtained for hot section components using these models. Future 
work is also discussed. 


VISCOPLASTIC MODELS 

Robinson* s models for isotropic materials and metal matrix composites are based 
on the concept of a flow potential. The flow and growth laws for the internal 
variables are derived from this flow potential. The material behavior is elastic for 
all the stress states within the potential, and it is viscoplastic for all the stress 
states outside the flow potential. The displacements (strains) are assumed to be 
small, and the total strain is assumed to be the sum of elastic and inelastic compo- 
nents. The models include two internal state variables to account for the kinematic 
and isotropic hardenings. 

A brief description of the constitutive equations of these two models and the 
values of constants utilized in the numerical computations are given below. 
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Isotropic Model 


The constitutive equations for the isotropic model developed by Robinson 
(ref. 1 ) are 

Flow law. - 
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In the above equations, a^j denotes the stress, the inelastic strain 

rate, Sj^j the Kronecker delta, and £ij the effective stress. Repeated 
subscripts imply summation. k is the drag stress and a^j is the back stress. 

Numerical values of constants . - The values of constants are for 2-l/4Cr-lMo 
steel . 


VI = 3 . 61xl0 7 
n = 4 


VI = >1 exp (©i ) 
m = 7.73 
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(5 = 1.5 R = 9.0x10 8 exp (0£ ) 

H = 1.37xl0“ 4 G 0 = 0.14 

and 

6 l = (23.80 - 2635) 

®2 - 40 000 (sir 

Here, 0 is the absolute temperature in degrees Kelvin. 



Metal Matrix Composite Model 

The constitutive equations for this model, taken from reference 2 
summarized as follows: 

Flow law. - 
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5 = 



C = 


4Cco 2 - 1) 
4co 2 - 1 


Here, denotes a unit vector along the preferred direction at a point of the 

material, and w and n denote the anisotropic parameters. 


Numerical values of constants . - In the absence of experimental data, the 
following values of the constants are assumed to carry out the numerical 
computations : 

p = 3.60x10^ n = 4 

m = 7.73 B = 0.75 


G 0 = 0.005 K t = 0.82 

R = 8.87xl0“ 8 w = 7 


n = 3 


IMPLEMENTATION IN MARC 

The user subroutine HYPELA in the MARC code provides the user with an efficient 
way of implementing and integrating the general nonlinear constitutive relationships 
of Robinson’s model. The details of this implementation may be found in reference 3. 

To avoid the numerical instabilities encountered due to the "stiff" and 
"discontinuous" nature of constitutive equations, the use of "smoothing functions" 

(or "spline functions") is made. The spline functions used in the computations are 
defined in references 1 and 3. 


RESULTS 

Numerical results using the Robinson’s models are depicted in figures 1 to 18. 
Figures 1 to 9 exhibit the results for the isotropic model, whereas figures 10 to 18 
show the results for the metal matrix composite model. 


Isotropic Model 

The cyclic thermal and mechanical loads used to generate the uniaxial hysteresis 
loops are shown in figure 1. Three different strain rates for isothermal, non- 
isothermal in-phase, and nonisothermal out-of-phase loadings are shown in figures 2 
to 4, respectively. Figure 5 shows excellent agreement between experiment and 
predictions, and verifies the correct finite element implementation of the model. 

The results for the multiaxial problem of a thick internally pressurized cylin- 
der are exhibited in figures 6 to 9. The analytical results shown in these figures 
were obtained using the expressions given in reference 3. The MARC and analytical 
results agree very well and establish the feasibility of performing a nonlinear 
structural analysis using the viscoplastic model. 
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Metal Matrix Composite Model 


The results of uniaxial calculations using Robinson's metal matrix composite 
model are depicted in figures 10 to 18. These figures show the hysteresis loops, 
creep curves, and relaxation curves generated with the finite element implementation 
of the model. The implementation clearly shows all the expected characteristics 
(e.g., strain-rate sensitivity, creep, relaxation, etc.) and the effect of preferred 
orientations (anisotropy) of the composites. The experimental work to determine the 
values of constants in the model for a metal matrix composite material is in progress 
at Lewis. 


CONCLUSIONS AND FUTURE WORK 

The viscoplastic models for the isotropic material and metal matrix composites 
developed by Robinson have successfully been implemented in MARC. The results 
obtained (in particular using the isotropic model) support the feasibility of 
performing nonlinear structural analyses using these models. These models will 
provide for more realistic and rational designs of gas turbine engine hot section 
components . 

The implementation of other viscoplastic models in MARC - such as (1) visco- 
plastic models for nickel-base superalloys, (2) viscoplastic models for single 
crystal alloys, and (3) viscoplastic models for metal matrix composites including 
damage - is in progress. The implementation of these models and their results, 
uniaxial and multiaxial predictions, will be reported subsequently. 
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CYCLIC THERMAL AND MECHANICAL LOADINGS 
USED TO GENERATE HYSTERESIS LOOPS 
BY FINITE ELEMENT CODE— MARC 
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Figure 1 


HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES 
USING ROBINSON’S MODEL— ISOTHERMAL CASE, 800 «F 
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Figure 2 
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HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES 
USING ROBINSON’S MODEL— NONISOTHERMAL CASE 

IN-PHASE THERMAL AND MECHANICAL LOADINGS 



Figure 3 


HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES 
USING ROBINSON’S MODEL— NONISOTHERMAL CASE 

OUT-OF-PHASE THERMAL AND MECHANICAL LOADINGS 
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Figure 4 




COMPARISON OF MARC AND EXPERIMENTAL 
HYSTERESIS LOOPS AT 1000 »F 



Figure 5 


STRESS DISTRIBUTION IN AN INTERNALLY 
PRESSURIZED THICK CYLINDER USING ROBINSON S MODEL 

(FEM-MARC SOLUTION) 
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Figure 6 
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COMPARISON OF MARC AND ANALYTICAL RESULTS 
FOR STRAIN DISTRIBUTION IN AN INTERNALLY 
PRESSURIZED THICK CYLINDER USING ROBINSON’S MODEL 
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Figure 7 


COMPARISON OF MARC AND ANALYTICAL RESULTS 
FOR STRESS DISTRIBUTION IN AN INTERNALLY 
PRESSURIZED THICK CYLINDER USING 
ROBINSON’S MODEL 
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Figure 8 
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STRAIN DISTRIBUTION IN AN INTERNALLY 
PRESSURIZED THICK CYLINDER 
(FEM-MARC SOLUTION) USING ROBINSON’S MODEL 
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Figure 9 


HYSTERESIS LOOPS AT DIFFERENT STRAIN 
RATES; ORIENTATION ANGLE = 0° 



Figure 10 
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HYSTERESIS LOOPS AT DIFFERENT STRAIN 
RATES; ORIENTATION ANGLE = 30° 



Figure 11 


HYSTERESIS LOOPS AT DIFFERENT STRAIN 
RATES; ORIENTATION ANGLE = 60° 
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Figure 12 
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HYSTERESIS LOOPS AT DIFFERENT STRAIN 
RATES; ORIENTATION ANGLE = 90° 



Figure 13 


CREEP CURVES FOR DIFFERENT FIBER 
ORIENTATION ANGLES 



Figure 14 
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RELAXATION CURVES FOR DIFFERENT FIBER 
ORIENTATION ANGLES 
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Figure 15 


HYSTERESIS LOOPS FOR DIFFERENT FIBER 
ORIENTATION ANGLES; STRAIN RATE = 0.001/min 
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Figure 16 
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HYSTERESIS LOOPS FOR DIFFERENT FIBER 
ORIENTATION ANGLES; STRAIN RATE = 0.01/MIN 
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Figure 17 


HYSTERESIS LOOPS FOR DIFFERENT FIBER 
ORIENTATION ANGLES; STRAIN RATE = 0.1/MIN 



Figure 18 
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